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1.  INTRODUCTION 


A.  GENERAL 

The  underwater  environment  is  the  most  challenging  region  of  the  battle  space  in 
which  to  communicate  effectively.  Sound  propagation  in  an  underwater  communication 
channel  becomes  a  very  difficult  task  due  to  numerous  constraints  and  limitations 
imposed  by  the  nature  of  the  medium  and  the  environment.  The  most  important 
characteristic  s  of  the  underwater  environment  that  put  limitations  on  the  performance  of 
underwater  communication  are  dispersion  caused  by  the  spatial  and  temperal  variability 
and  multipath  propagation.  Multipath  propagation  produces  intersymbol  interference 
(ISI)  and  severe  degradation  of  the  bit  error  performance  in  digital  communication.  The 
specific  investigation  in  this  thesis  is  the  demodulation  of  Binary  Phase- Shift  Keying 
(BPSK)  in  three  different  shallow  water  ocean  acoustic  channels. 

B.  OBJECTIVE 

The  objective  of  this  thesis  was  to  study  the  possibility  of  using  the  time-reversed 
ocean  impulse  response  prior  to  the  demodulation  of  a  BPSK  signal  to  mitigate  the 
degradation  of  the  bit  error  performance  in  multipath  propagation.  The  thesis  evaluates 
the  bit  error  performance  of  the  mitigated  signal  compared  to  the  bit  error  performance  of 
the  received  signal  distorted  by  the  intersymbol  interference  in  three  different  ocean 
environments  in  the  absence  of  noise. 

C.  ORGANIZATION 

The  remainder  of  this  thesis  is  organized  as  follows.  Chapter  II  gives  the  general 
theory  related  to  BPSK  binary  communication.  This  describes  the  definition,  waveform 
and  constellation  of  the  BPSK  signal,  the  theoretical  convergence  of  the  bit  error 
probability,  and  two  types  of  the  bit  error  degradation  factors.  Chapter  III  describes  the 
experiments  and  simulation  procedure  used  in  this  thesis.  This  includes  a  time- domain 
and  frequency- domain  BPSK  representation.  Power  Spectral  Density  (PSD)  for  BPSK,  a 
description  of  the  Monterey- Miami  Parabolic  Equation  (MMPE)  model,  the  ocean 
impulse  response  extraction  from  MMPE  model,  BPSK  demodulation  and  multipath 
mitigation.  Chapter  IV  provides  the  main  results  of  the  simulation.  Included  in  the  results 
are  the  evaluation  for  the  bit  error  performance  of  the  BPSK  signal,  bit  error  degradation 
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and  multipath  mitigation  for  BPSK  signal  in  three  different  shallow  water  acoustic 
channels. 
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n.  BINARY  PHASE-SHIFT  KEYING 


This  chapter  describes  the  general  theory  related  to  binary  phase- shift  keying 
(BPSK).  First,  the  definition,  waveform  and  its  constellation  for  BPSK  modulation  are 
explained.  Next,  as  an  important  criterion  of  performance  in  digital  communication,  the 
probability  of  bit  error  is  defined  and  the  theoretical  convergence  of  the  bit  error 
probability  is  described.  Finally,  two  types  of  the  bit  error  degradation  factors  are 
presented. 

A.  BINARY  PHAS&SHIFT  KEYING 

In  binary  phase-shift  keying  (BPSK),  the  phase  of  a  constant- amplitude  carrier 
signal  is  switched  between  two  values  according  to  the  two  possible  “messages”  and 
corresponding  to  binary  1  and  0,  respectively.  Normally,  the  two  phases  are  separated 
by  180  degrees.  If  the  sinusoidal  carrier  has  amplitude  A^ ,  then  the  average  energy  per  bit 
1  2 

is  E/,  =  —  A^  where  T,,  is  the  bit  time  duration,  and  the  transmitted  BPSK  signal  can 
be  represented  as: 

s,(t)=JV^cos[27l//+0,(f)]  .  ,  ‘  (2.1) 

V  T*  «  =  1,2 


where  the  phase  term,  0,  (0 ,  will  have  two  discrete  values,  given  by  0j=O  and  02  =7t  . 
In  other  words,  in  BPSK  modulation,  the  modulating  data  signal  shifts  the  phase  of  the 
waveform  5,.(t)  to  one  of  two  states,  either  zero  or  71(180°).  Instead  of  using  a 

rectangular  shaped  pulse  as  shown  in  Figure  2.1,  we  can  use  non- rectangular  shaped 
pulses  by  applying  banning,  hamming,  and  Gaussian  windows.  However,  we  shall  adhere 
to  the  simplest  case  of  a  rectangular  pulse  in  this  thesis. 


Figure  2.1.  BPSK  Signal  in  the  Time  Domain. 
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Figure  2.2.  BPSK  Signal  in  Frequency  Domain. 

Figure  2.1  shows  a  typical  BPSK  waveform  in  the  time  domain  with  its  abrupt 
phase  changes  at  the  symbol  transitions.  The  sketch  in  Figure  2.2  shows  the  spectrum  of 
the  BPSK  waveform  in  frequency  domain.  For  this  signal  set,  the  single  waveform 


=  4^  cos(271  fj)  0  <  t  < 


(2.2) 


can  be  defined.  Using  this  signal,  the  BPSK  signal  set  can  be  represented  as 


^BPSK  ~  {  4^* 


(2.3) 
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B.  SIGNAL  CONSTELLATION  FOR  BPSK 


Vector 


M  =  2 


S2 


Si 


Figure  2.3.  Signal  Constellation  for  BPSK. 


A  BPSK  signal  waveform,  as  shown  in  Figure  2.3,  can  be  geometrically 
represented  as  vectors  or  phasors  on  a  polar  plot;  the  vector  length  corresponds  to  the 
signal  amplitude,  and  the  vector  direction  corresponds  to  the  phase.  Such  a  representation 
is  referred  to  as  the  symbol  constellation  and  provides  a  graphical  representation  of  the 
complex  envelope  of  a  BPSK  signal  for  the  two  possible  symbols.  The  taxis  of  the 
constellation  diagram  represents  the  “in-phase”  component  of  the  complex  envelope,  and 
the  Q-axis  represents  the  “quadrature”  component  of  the  complex  envelope.  The  distance 
between  the  signals  on  a  constellation  diagram  measures  separation  of  the  modulation 
waveforms  and  determines  how  well  a  receiver  can  differentiate  between  all  possible 
symbols  when  random  noise  is  present.  The  larger  the  signal  distance,  the  better  the 
chance  of  correct  symbol  detection. 

C.  PROBABILITY  OF  BIT  ERROR  FOR  BPSK 

An  important  measure  of  performance  used  for  digital  modulation  is  the 
probability  of  symbol  error,  Pg .  It  is  often  convenient  to  specify  system  performance  by 
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the  probability  of  bit  error  Pg  ,  even  when  decisions  are  made  on  the  basis  of  symbols 
rather  than  bits.  The  relationship  between  Pg  and  Pg  for  orthogonal  signaling  is: 


P 


B 


Pe 

log2M 


(2.4) 


where  M  is  the  number  of  symbols. 

For  BPSK  modulation  (M=2),  the  symbol  error  probability  is  equal  to  the  bit  error 
probability.  When  the  signals  are  assumed  equally  likely  and  signal  (j  =  l,2)  is 

transmitted,  the  received  signal,  q{t) ,  is  equal  to  s.  (t)  +  n(t) ,  where  n(t)  is  modeled  as 
additive  white  Gaussian  noise  (AWGN).  The  antipodal  signals  (signals  of  equal 
amplitude  and  opposite  polarity)  s^(t)  and  are  [cf.  Eq.  (2.3)]: 


i(0  =  V^vVi(0  0<t<T^ 

(2.5) 

(0  =  0  <  t  <  T^ 

(2.6) 

The  detection  and  demodulation  of  the  received  signal  involves  a  correlation  with 
each  of  the  antipodal  signals,  as  explained  in  Chapter  3.  At  the  end  of  each  symbol 
duration  T^ ,  the  output  of  the  correlator  yields  a  sample  z(T^),  called  the  detection 
statistic  given  by 


z(Tj  =  v,(T.)  +  n„(Tj  i  =  l,2  (2.7) 

where  v,  ( T^ )  is  the  desired  signal  component  and  ( T^ )  is  the  noise  component.  The 
noise  component  is  a  zero  mean  Gaussian  random  variable,  and  thus  z  (T^ )  is  a  Gaussian 
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random  variable  with  a  mean  of  either  Vj  or  Vj  depending  on  whether  (?)  or  (0 
sent. 

The  decision  stage  of  the  detector  will  choose  the  with  the  largest  correlator 
output  Zi  (?) ,  or  in  this  case  of  equal-energy  antipodal  signals,  the  detector,  using  the 
decision  rules,  decides: 


5i(?)  if  z(Tt)>lo 

(2.8) 

S2  (?)  otherwise 

(2.9) 

where  denotes  the  decision  threshold  (equal  to  0  for  equally  probable  antipodal 
signals)  and  z(T^)  is  the  correlator  output  at  time  .  Two  types  of  detection  errors  can 
be  made.  The  first  type  of  error  takes  place  if  jj(?)  is  transmitted  but  the  detector 
measures  a  negative  value  for  z(T^)  and  decides  (incorrectly)  that  signal  ^2(0  was  sent. 
The  second  type  of  error  takes  place  if  signal  ^2(0  is  transmitted  but  the  detector 
measures  a  positive  value  for  z(T^)  and  decides  (again  incorrectly)  that  signal  5j(?)  was 
sent.  Therefore,  the  probability  of  bit  error,  Pg  ,  is  the  sum: 

PB  =  ^(^2ki)^(^i)  +  ^(^ik2)P(^2)  (2-10) 

where  )  are  the  conditional  probabilities  and  P(5,)  are  the  prior  probability  of  the 

symbols.  For  the  case  when  the  symbols  are  equally  likely  (which  is  mostly  the  case): 

P(s,)=P(i,)  =  i  (2.11) 

The  expression  for  the  bit  error  probability  then  becomes: 
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(2.12) 


PB=^^’(j2kl)+^^’(^LK) 

The  conditional  probabilities  P(52|'5'i)  and  are  found  by  integrating  the 

conditional  probability  density  function  (pdf)  of  the  output  of  the  correlator  which  is 
depicted  in  Figure  2.4 


Figure  2.4  Conditional  Probability  Density  Functions:  p[^z\s^), 

The  conditional  probability  density  functions  (pdf’s)  of  the  signal  with  an  additive 
white  Gaussian  noise  (AWGN)  are: 

1  r  r  -  t” 

p[z\s,)  =  —i= — exp  -  i  =  l,2  (2-13) 
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i.e.,  the  conditional  pdfs  p{z\s.)  {i  =  l,2)  are  Gaussian  with  mean  value  v,,  and 

2  2 

variance  a  ^  where  c  ^  corresponds  to  the  noise  variance  at  the  output  of  the  correlator. 
Because  of  the  symmetry  of  p{^z\s.) ,  the  bit  error  probability  of  Eq.  ( 2. 12)  reduces  to: 

(2-14) 

Thus,  the  probability  of  a  BPSK  bit  error  Pg  is  numerically  equal  to  the  area 
under  the  “tail”  of  either  pdf  p(z|5j)  or  p  (z|5’2)  that  falls  on  the  incorrect  side  of  the 
threshold.  We  can  therefore  compute  Pg  by  either  integrating  p(z|5j)  between  the  limits 
-oo  and  Yg  or  by  integrating  p(z|52)  between  the  limits  and  0°,  where 
Yo  =  ( Vj  +  V2  )/2  is  the  optimum  threshold.  Hence, 


Pb  = 

7o 


(2.15) 


Using  Eq.  (2.13) ,  the  probability  of  bit  error  for  BPSK  is: 


r  ^ 

1 

^  z  -  V2  ^ 

2 

~2 

1  ^0  J 

dz 


(2.16) 


If  we  introduce  u  - - -  and  o^du  =dz,  the  integral  simplifies  to 


I  .R 


g  j  p. — .-  exp 


V  ^  J 


du  =  Q 


(  ..  3 


V  y 


(2.17) 
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where  Q{-)  is  the  complementary  error  function  defined  as: 


Q(x)  = 


.2  \ 


du 


(2.18) 


For  equal  energy  antipodal  signaling,  such  as  the  BPSK  format  in  Eq.  (2.3),  the 

receiver  output  signal  components  are  Vj  =  when  5j  (t)  is  sent  and  V2  =  when 

2 

52(0  is  sent.  For  AWGN,  the  noise  variance  of  the  correlator  output  can  be  written 
as  /2,  where  /2  is  two-sided  power  spectral  density  of  the  noise.  Thus,  we  can 
rewrite  probability  of  bit  error  as  follows: 


Pb=  ^ 


-J2E,/N„ 


exp 


.2  3 


\iu  =  Q 


(2.19) 


D.  BIT  ERROR  PROBABILITY  CONVERGENCE 

In  digital  communications,  we  more  often  use  E^/A^  ,  a  normalized  version  of 

average  signal  power  to  average  noise  power  ratio  (S/N  or  SNR ),  as  a  figure  of  merit. 
E^  is  bit  energy  and  can  be  described  as  signal  power  S  times  the  bit  time  is 

noise  power  spectral  density,  and  can  be  described  as  noise  power  N  divided  by 
bandwidth  W .  Since  the  bit  time  and  bit  rate  are  reciprocal,  we  can  replace  with 
l/R^  and  write 


Pfc  _  (220) 

N/W  N/W 

If  we  use  R  instead  of  R^  to  represent  bits/sec,  and  we  rewrite  Equation  (2.20) 
to  emphasize  that  E^/A^  is  just  a  version  of  S/N  normalized  by  bandwidth  and  bit  rate, 
we  obtain: 
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(2.21) 


n^~n[r^ 

One  of  the  most  important  metrics  of  performance  in  a  digital  communication 
system  terms  is  a  plot  of  the  bit  error  probability  Pg  versus  .  Figure  2.5  shows  the 

“waterfall- like”  shape  of  most  such  curves. 


Pb 


For  E^/A^q  >  v:Q,Pg  <Pq.  The  dimensionless  ratio  E^/A^^  is  a  standard  quality 
measure  for  digital  communications  system  performance.  Therefore,  the  smaller  the 
required  E^  /  A^^  ,  the  more  efficient  is  the  detection  process  for  a  given  probability  of 
error. 
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To  optimize  (minimize)  Pg  in  the  context  of  an  AWGN  channel  and  the  receiver, 
we  need  to  select  the  optimum  receiving  filter  and  the  optimum  decision  threshold.  For 
the  binary  case,  the  optimum  decision  threshold  is  Yo  =  (^i  +  ^2)/2  was  shown  in 

Equation  (2.17)  that  this  threshold  results  in  Pg  =  e[(vj-vJ/2aJ.  Next,  for 
minimizing  Pg ,  it  is  necessary  to  maximize  the  argument  of  Q{x)  .  Thus,  we  need  to 
maximize  (vj  -  Vj  )/2cJo  or  equivalently,  maximize 


(v.  -V, )  ,  . 

=  ^ ^  (2.22) 

where  (vj  -Vj)  is  the  difference  of  the  desired  signal  components  at  the  filter  output  at 

time  t  =  T^ ,  and  the  square  of  this  difference  signal  is  the  instantaneous  power  of  the 

2  2 

difference  signal.  The  5AY?  can  also  be  maximized  by  minimizing  ;  however,  is 
usually  not  under  our  control. 

E.  TWO  TYPES  OF  BIT  ERROR  PERFORMANCE  DEGRADATION 

The  effects  of  bit  error  performance  degradation  can  be  partitioned  into  two 
categories.  The  first  is  due  to  a  decrease  in  received  signal  power  or  an  increase  in  noise 
or  interference  power,  giving  rise  to  a  loss  in  signal- to- noise  ratio  or  E^/A^  .  The  second 
is  due  to  signal  distortion,  such  as  might  be  caused  by  intersymbol  interference  (ISI) 
where  the  tail  of  a  pulse  can  “smear”  into  adjacent  symbol  intervals,  thereby  interfering 
with  the  detection  process  and  degrading  the  error  performance.  Even  in  the  absence  of 
noise,  the  effects  of  filtering  and  channel- induced  distortion  lead  to  ISI.  Eigure  2.6a 
shows  the  effect  of  a  loss  in  E^/A^;  the  solid- line  corresponds  to  the  theoretical  Pg 

versus  E^/A^  curve  and  the  dashed- line  corresponds  to  a  degradation  effect  brought 
about  by  a  loss  in  E^/A^  .  Eigure  2.6b  shows  an  irreducible  Pg  caused  by  distortion 
(ISI).  The  solid- line  corresponds  to  the  theoretical  Pg  versus  E^/A^  curve  and  the 
dashed-  line  corresponds  to  a  degradation  effect  brought  about  by  ISI  instead  of  a  simple 
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loss  in  signal- to- noise  ratio.  The  bit  error  rate  in  this  case  cannot  further  be  reduced  by 
increasing  the  signal  power. 


Figure  2.6.  (a)  Loss  in  .  (b)  Irreducible  Pg  Caused  by  Distortion 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 
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m.  EXPERIMENTAL  DESCRIPTION 


In  this  chapter,  we  describe  the  experiments  and  simulation  procedure  used  in  this 
thesis.  First,  as  a  part  of  BPSK  signal  generation,  time-domain  and  frequency- domain 
BPSK  representations  and  Power  Spectral  Density  (PSD)  for  BPSK  are  analytically 
derived.  Second,  we  describe  the  Monterey- Miami  Parabolic  Equation  (MMPE)  acoustic 
simulation  model  and  how  to  generate  the  ocean  responses  in  the  time- domain  and 
frequency- domain  from  the  MMPE  model.  Third,  we  describe  the  ocean  environment  for 
our  simulations  and  the  input  data  for  the  MMPE  model.  Einally,  we  present  the  results  of 
BPSK  demodulation  and  multi- path  mitigation. 

A.  BPSK  SIGNAL  GENERATION 

In  Chapter  II,  we  described  the  properties  of  BPSK.  In  this  section,  time- domain 
and  frequency- domain  representations  for  the  BPSK  waveform  are  described  and  related 
to  the  ocean  impulse  response,  and  the  power  spectral  density  of  the  received  signal. 
Eigure  3.1  shows  a  diagram  for  BPSK  signal  generation.  The  quantity  b{t)  is  a  random 

binary  signal  at  baseband.  The  message  data,  {a„}  represent  a  set  of  random  variables 


with  =  ±1  and  probabilities  P(a„  =  l)=  P(a„  =  -l)=  The  function  g{t-nT^)  is  a 

shaped  pulse  in  general;  however  we  will  take  it  to  be  a  rectangular  pulse  with  bit  time 
duration,  T^.  The  signal  ^(t)  is  a  sinusoidal  carrier  with  center  frequency  and 


amplitude  .  The  output  s(t)  is  the  modulated  BPSK  signal. 


^(0=  £  a^g{t-nTj 

n=-o°  ^ 

^  ►  s(t)=b(t)-<^(t) 

i. 

Random  binary  signal 

BPSK  signal 

4)(0=Acos(27r/,f) 

Sinusoidal  carrier 

Eigure  3.1.  BPSK  Signal  Generation. 


15 


The  Fourier  transform  for  a  rectangular  pulse  with  amplitude  one  and  bit  duration 


is  given  by: 


G(/)  = 


sin(7t/TJ 


(3.1) 


The  Power  Spectral  Density  (PSD)  for  a  random  binary  signal,  bit),  will  be 
evaluated  by  first  truncating  bit)  as  follows: 


-Tj2<t<Tj2 

n=-N 


(3.2) 


where  bj  it)  is  the  truncated  signal  and  T^  /2  =  (A^  +1/2)  T^-  The  Fourier  transform  of  the 
truncated  signal,  bjit),  is  then. 


Brif)=FT[brit)h  £a„Fr[g(t-nTj]=  Y,a„Gif)e-^^^^ 

n=-N  n——N 

=  (3.3) 

n=-N 


where  the  notation  FT[]  represents  the  Fourier  transform  operation  and  G{f)  is  the 
resulting  Fourier  transform  of  g  (t) . 

The  PSD  for  the  random  signal  bit)  can  be  defined  by  [Ref  3] 


To^” 


16 


where  Bj  (/ )  represents  the  Fourier  transform  of  the  corresponding  truncated  signal  and 
the  overbar  represents  statistical  expectation.  Upon  substituting  (3.3)  into  (3.4),  the  PSD 
of  b{t)  becomes 


Art 


-j(an\ 


n=-N 


/  1  w  3 

=  |G(/f  lim  — Vl  =|g(/)|  lim 


2N  +  1 
'  T 


where  we  have  assumed  that  the  are  statistically  independent.  Now,  using 
=  2 ( + 1/2) ,  we  obtain 


2N  +  1 

(2M  +  \X 


T, 


j-iG(/)r=T, 


sin(7t/rj 


(35) 


where  Eq.  (3.l)  has  been  used  for  G(/) .  This  is  the  PSD  of  the  signal  at  “baseband.’ 
For  the  modulated  BPSK  signal  we  have 


s(t)  =b(t)A^  cos{2nfj) 


(3.6) 


The  corresponding  PSD  is  given  by 

.2  .2 

^.(/)=^n(/  +  /J+^n(/-/J  (3.7) 

4  4 

(see  Ref  3),  where  is  the  PSD  of  the  baseband  signal  b(t) .  Upon  substituting  (3.5) 
into  (3.7)  we  obtain 
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(3.8) 


Ps(f)  = 


4 


Pif  +  +\G{f  - 


sin7l(/  +  /jT, 

2 

sin7l(/  -/JT, 

2' 

- < 

4 

Af+lX  J 

+ 

J 

B.  MMPE  MODEL  DESCRIPTION 

The  Parabolic  Equation  (PE)  method  was  introduced  into  underwater  acoustics  in 
the  early  1970’ s  by  Tapper!  [Ref  8].  The  Monterey- Miami  Parabolic  Equation  (MMPE) 
model  [Ref  9]  is  a  numerical  model  to  solve  acoustic  wave  propagation  problems  in  the 
ocean  ising  the  PE  method.  The  MMPE  Model  is  used  to  compute  the  response  at  the 
receiver  location  of  a  signal  transmitted  from  a  source  and  traveling  through  the  ocean. 
Eigure  3.2  shows  the  three  stages  of  using  this  model.  In  the  first  stage,  we  characterize 
the  input  parameters;  in  the  second,  we  run  the  model;  and  in  the  last  stage,  we  extract  the 
frequency  response  and  ocean  impulse  response  from  the  outputs  of  the  MMPE  model. 


Stage  1: 

^ - ^ 

Stage  3: 

Characterize  the 

Stage  2: 

Extract  the  Ocean 

Input  Parameters 

—J,  Run  MMPE  ) _ ^ 

Erequency  Response  and 

for  MMPE 

model 

Ocean  Impulse  Response: 

model 

^ - ^ 

H(f)  and  h(t) 

Eigure  3.2.  Three  Stages  to  Extract  Ocean  Response  from  MMPE  Model 

1.  A  Brief  Description  of  the  MMPE  Model 

In  order  to  derive  the  frequency- domain  and  time- domain  response  of  the  ocean 
underwater  communication  channel,  the  general  theory  behind  the  parabolic  equation 
model  will  be  introduced.  Its  implementation  -  the  split- step  Eourier  algorithm  for  solving 
the  PE  (PE/SSE)  method  [Ref  10]  is  also  briefly  described. 

In  ocean  communication,  all  signals  such  as  p(t)  represent  acoustic  pressure  as  a 
function  of  time,  and  are  typically  denoted  by  p{x,t)  where  v  represents  the  ocean 
coordinates  at  which  the  pressure  signal  is  measured.  The  quantity  p[x,t)  is  often 
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referred  to  as  a  space- time  signal,  since  it  is  a  function  of  the  spatial  parameters  x  and 

time  t .  When  measured  at  a  fixed  (receiver)  location  (say  x  =  x\)  it  becomes  a  function 
of  time  only.  It  is  transformed  into  an  electrical  signal  by  an  acoustic  transducer  or 
hydrophone. 

The  inhomogeneous  wave  (Helmholtz)  equation  for  the  acoustic  pressure  p(x,t) 
in  a  medium  with  sound  speed  and  density  p  (xj  can  be  expressed  as 


where  represents  the  sound  source  distribution  as  a  function  of  both  position  x 

and  time  t .  In  our  application  (as  in  most)  there  is  a  single  source  at  a  particular  location 
x=Xo.  Hence  5(x,i)  is  formally  a  spatial  “impulse”  at  Xo  with  time  dependence 
s{t)  corresponding  to  the  transmitted  waveform. 

The  parabolic  equation  model  is  based  on  an  approximation  of  the  Helmholtz 
wave  equation  in  a  cylindrical  coordinate  system.  Because  of  the  ocean’s  relative 
shallowness  compared  to  horizontal  propagation  distance  for  the  majority  of 
environments,  it  is  well  suited  for  a  description  in  cylindrical  coordinates. 

The  signals  of  interest  in  underwater  communications  such  as  BPSK  have 
complicated  time  dependence;  thus  analytical  solution  to  the  wave  equation  (3.9)  would 
be  impossible.  However  the  equation  is  linear,  and  a  solution  can  be  computed  relatively 
easily  for  sources  that  have  time  dependence  of  the  form  e  for  a  single  frequency  /  . 

Thus  for  a  source  with  a  time  dependence  of  the  form  ,  the  output  at  a  point  x  =  xi 

is  dependent  on  /  and  can  be  written  as  H  ■  (/)e^^"^'.  This  is  referred  to  by  people 

working  in  ocean  acoustics  as  a  ‘time  harmonic”  solution  to  the  wave  equation  The 
quantity  H  -  (f)  is  the  ocean  frequency  reponse  and  is  inverse  Fourier  transform  is  the 
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ocean  impulse  response  (Green’s  function).  The  time  harmonic  solution  forms  the  basis 
of  computing  the  response  due  to  a  source  with  arbitrary  time  dependence  by  using 
Fourier  analysis/synthesis  methods. 

From  this  point  forward  the  analysis  will  be  carried  out  in  discrete  time.  This  is 
appropriate  since  numerical  techniques  (i.e.,  the  DFT)  are  used  to  perform  the 
computations  in  both  the  time  and  frequency  domains.  We  focus  first  on  the  analysis  of 
low-pass  signals  at  baseband  —W/2<f<Wj2  and  later  extend  the  discussion  to  the 
case  of  bandpass  signals  such  as  BPSK  centered  around  some  higher  carrier  frequency.  In 
particular,  let  us  sample  a  signal  s(t)  at  times  t  =  0,  T ,  2T and  represent  the 
resulting  discrete-time  signal  s\p\  =  s{nT)  in  terms  of  its  Discrete  Fourier  Transform 
(DFT): 


s 


.2%  kn 


•  k=0 


(3.10) 


Here  .sfl;]  is  the  DFT  eoefficient  representing  the  frequency  component  of  the  signal  at 

frequency  f  =  klLT : 


L-l 


.2%kn 


n=0 


(3.11) 


and  L  is  the  number  of  time  and  frequency  samples  used  in  the  transform.  The  number  of 
samples  L  and  time  interval  T  must  of  course  be  chosen  to  cover  the  time  duration  and 
bandwith  of  the  signal. 

With  the  representation  Eq.  (3.10)  for  the  source  waveform,  the  output  waveform 
at  a  point  xi  can  be  computed  from 
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(3.12) 


'  k=0 


\  ^  J 


where  the  ocean  frequency  response  has  been  evaluated  at  the  frequency  /  =  27r^/L  and 
multiplied  by  the  Fourier  coefficient  5[^]  for  the  source. 

The  ocean  frequency  response  terms  that  need  to  be  computed  by  MMPE  are 


H- 


2Kk\ 


k^O,  1,2,  ...,  --1 


(3.13) 


which  represent  frequencies  up  to  the  Nyquist  “fold-over”  frequency  /^/2  =  ^^.  The 

other  terms  in  the  transform  represent  the  “negative”  frequencies  and  may  be  computed 
from 


fink  ^ 

=  H- 

'2%{L-ky 

[  ^  ) 
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(3.14) 


Both  the  forward  DFT  in  Eq.  (3.1l)  needed  to  compute  the  frequency  coefficients  5[^] 
of  the  source  signal  and  the  inverse  DET  represented  in  Eq.  (3.12)  are  performed  using 
an  EET  program. 

To  compute  H  ,{f)  we  proceed  as  follows.  Assuming  a  time  harmonic  solution, 
the  Helmholtz  equation  takes  the  form 
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(3.15) 


1  3  r  3p(r  ^  tp  ,/)  ^ 
rdr[^  df  ^ 


1  3V(r,z,(p,/)  ^ 

acp^ 


d^pjr  z  P  /) 
az" 


’9)P(  r,ZfS?,f)=  -471  PqS  (r  -  rj, 


where 


p(r  .z  IP  ,/)  =  /j(r,zip>? 


(3.16) 


The  reference  wavenumber  is  related  to  a  reference  sound  speed,  Cq  ,  by 


k 


0 


271/ 


and  the  acoustic  index  of  refraction  is  defined  by 


1  = 


c 


0 


c(r,z 


(3.17) 


(3.18) 


Note  that  in  this  derivation  the  density  variances  are  neglected.  By  defining  the  effective 
index  of  refraction  which  contains  the  appropriate  addition  terms,  the  influence  of  the 
density  differences  at  the  water-bottom  interface  can  be  included  [Ref  8]. 

By  assuming  that  the  ocean  acts  as  a  waveguide  with  a  cylindrical  coordinate 
system,  acoustic  energy  is  mainly  propagated  outward  from  a  source  in  the  horizontal 
direction.  Therefore,  the  pressure  field  can  be  approximated  by 


p(r,z,(p,/)  =  'T(r,z,(p,/)H 


(3.19) 
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where  'P(r, z,(p,/)  is  a  slowly  varying  envelope  function,  or  PE  field  function  and 
H  is  the  zero-th  order  Hankel  function  of  the  first  kind  of  the  outgoing  acoustic 

wave.  Taking  advantage  of  the  far- field  (kor»l)  asymptotic  approximation  of  the 
Hankel  function,  Eq.  (3.19)  can  be  rewritten  as 


p(r,z,(p,/) 


4~r 


(3.20) 


normalized  such  that  at  the  reference  range  r  =  R^,  =  Pq.  Substituting  Eq.  (3.20)  into 

Eq.  (3.15)  and  dropping  the  source  term  on  the  right  hand  side  gives 


- 1 — 1“  - 1 r - - r — h 

dr^  dr  r^  dtp  dz^ 


(1^-1)  +  ^ 

V  )  4^2 


(3.21) 


Neglecting  the  azimuthal  coupling  and  the  far- field  terms,  and  dropping  the  first  term  due 
to  the  slow  modulation  of  the  envelope  function,  Eq.  (3.2l)  can  be  rewritten  as 


kg  dr 


1  d^T^ 

2kQ  dz^ 


(3.22) 


Now,  by  defining  the  operators 


1 

T  = - 

2h 


^  dz^ 


(3.23) 


and 
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(3.24) 


U 


op 


Eq.  (3.22)  can  be  rewritten  as 


k„  dr 


=  (r  +u  )'P  . 

\  op  op  J 


(3.25) 


Equation  (3.25)  is  known  as  the  “standard”  parabolic  equation  (SPE)  [Ref  8],  with 

aecurate  solutions  limited  to  a  half  width  of  ±15°  for  the  propagation  angle  and 
represents  the  eomplete  description  of  the  outgoing  acoustic  energy  propagation  in  the 
waveguide.  In  order  to  extend  this  limit  to  ±  40° ,  a  higher  order  wide-angle  parabolic 
equation  (WAPE)  approximation  [Ref  11]  is  used.  Its  operators  are  defined  by 


T'wAPE  ~ 


J_i! 


(3.26) 


and 


^WAPE  (^  1)  • 


(3.27) 


The  MMPE  uses  the  split- step  Eourier  (SSE)  method  in  order  to  solve  the 
parabolic  equation  numerieally.  This  algorithm  integrates  the  solution  in  range  by 
applying  the  and  the  Uy^/^pp-  operators  in  the  wavenumber  (k^  )-domain  and  the 
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spatial  (z  )-domain,  respectively,  where  each  operator  is  a  scalar  multiplier.  In  the  k^  - 


A 

space,  the  wide  angle  T  wape  operator  to  simplify  calculation  is  defined  as 


A 

T  WAPE 


(3.28) 


The  field  function  at  range  r  +  Ar  is  expressed  as 


'T(r  +  Ar, 


•  f\- 


Ar, 
-jkn — t 


i^o~^  wape[<' 


(3.29) 


Since  the  output  data  of  the  MMPE  model  is  the  field  function  in  the  form  of 
equation  (3.29),  we  substitute  Eq.  (3.29)  into  Eq.  (3.20).  The  ocean  frequency  response 
(with  respect  to  pressure)  at  range  r  and  depth  z  is  then  given  by 


(3.30) 


As  described  earlier,  for  broadband  signals,  the  output  results  require  knowing 
^{r,z)(f)  frequencies  in  the  bandwidth  of  the  signal  [see  Eq.  (3.12)].  Thus  the 

MMPE  model  must  be  run  for  all  frequencies  within  the  band.  The  model  produces 
frequency  domain  results  which  are  written  to  an  output  file.  A  post  processing  program 
performs  the  synthesis  specified  in  Eq.  (3.12). 

2.  Signal  Representation  in  MMPE 

The  time  intervals  encountered  in  ocean  acoustics  can  be  large.  In  particular  the 
time  due  to  propagation  delay  from  source  to  receiver  driving  which  no  sound  is  heard  at 
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the  receiver  can  be  very  large  compared  to  the  duration  of  the  signal  when  it  finally 
arrives  (see  Figure  3.3) 


Received 

Signal 


Figure  3.3.  Pulse  Arriving  after  Propagation  Delay  x  . 

To  obtain  the  most  accurate  results,  the  analysis  time  needs  to  be  focused  in  a  small  time 
window  surrounding  the  received  signal  as  shown  in  Figure  3.3. 
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The  signals  of  interest  are  also  often  restricted  in  frequency  as  well  as  time.  For 
example,  a  BPSK  transmission  has  a  frequency  spectrum  similar  to  that  shown  in  Figure 

3.4. 


Figure  3.4.  Spectrum  of  Bandpass  Signal 


Such  signals  are  referred  to  as  bandpass  signals;  the  spectral  power  is  assumed  negligible 
outside  of  an  interval  ±W f 2  around  the  center  frequency  or  carrier  frequency  .  Again, 
for  accuracy  and  reduced  computation  in  processing,  it  is  most  appropriate  to  focus  in  a 
frequency  window  (frequency  band)  where  the  spectrum  is  non- zero. 

MMPE  deals  with  the  signal  processing  in  an  efficient  way  by  permitting  the  user 
to  specify  a  time  window  and  frequency  band  for  the  processing.  The  formulation  for 
such  processing  is  described  below. 

a.  Time  Windowing 

The  time  windowing  in  MMPE  is  relatively  simple  to  explain.  The 

i2jt  /— 

frequency  response  of  Eq.  (3.30)  contains  a  term  e^^''=e  which 
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corresponds  to  the  gross  propagation  delay  x  =  r/c^  from  source  to  receiver  in  the  ocean. 
MMPE  drops  this  term,  so  that  the  modified  frequency  response  is  given  by 


(3.31) 


The  time  domain  signals  resulting  from  this  modified  frequency  response  correspond  to  a 
reduced  time  scale  t  defined  as 


'  r 

t  -t-x-t - 


(3.32) 


and  referred  to  as  reduced  time .  MMPE  uses  this  reduced  time  for  its  output  results. 
Absolute  time  (relative  to  the  source)  can  be  obtained  by  adding  back  in  the  gross 
propagation  delay  x  =  r/c^ . 

Without  the  modification  to  the  frequency  response  specified  in  Eq. 
(3.31) ,  some  time  response  would  be  produced  within  the  analysis  window  of  the  DET. 
However,  this  response  would  appear  as  a  circularly  shifted  version  of  the  true  correct 
response  and  the  time  scale  would  not  be  correct.  This  effect  is  known  as  ‘time  aliasing.” 
b.  Processing  of  Bandpass  Signals 

A  continuous -time  real  bandpass  signal  such  as  that  depicted  in  Eigure  3.4 
can  be  written  in  the  general  form 


(0  =  (0  cos27t  ft  -  s^  (t)  sin2jt  ft 


28 


Here  Sp{t)  and  s^{t)  are  known  as  the  “in-phase”  and  “quadrature”  modulation 
components.  The  bandpass  signal  can  be  equivalently  expressed  as 


s^{t)=Re 


(3.33) 


where 


^(0  =  ^;,  (0+i\(0 


is  known  as  the  complex  envelope  .  The  complex  envelope  is  a  complex- valued  low-pass 
signal  with  spectrum  as  depicted  in  Figure  3.5. 


Figure  3.5  Spectrum  of  a  Complex  Valued  Lowpass  Signal  Corresponding  to 

Bandpass  Signal  of  Figure  3.4. 
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The  complex  envelope  may  be  obtained  from  the  original  bandpass  signal  by  discarding 
the  negative- frequency  part  of  the  spectrum,  moving  the  positive- frequency  part  of  the 
spectrum  from  a  center  point  of  /  =  /^  to  /  =  0  (“baseband”)  and  scaling  by  a  factor  of 

2. 

The  advantage  of  using  complex  envelope  representation  for  the  bandpass 
signal  is  that  the  processing  can  be  focused  in  the  bandpass  region  by  processing  the 
complex  envelope  at  baseband.  Let  us  denote  the  sampled  version  of  the  complex 

envelope  as 


s[n]  =  s[nT)  =  Sp  (nT)  +  js^  [nT] 


(The  notation  ^[n]  was  previously  used  to  denote  a  sampled  real  band- limited  baseband 
signal;  the  only  difference  now  is  that  ^[n]  is  complex.)  The  complex  baseband  signal 
can  be  represented  as  in  Eq.  (3.10)  where  S[k]  are  the  frequency  components  at 
baseband  defined  by  Eq.  (3.11) .  However  a  frequency  in  the  DET  of  27r^/L  represents 
an  actual  frequency  in  the  ocean  of  InkfL+f^.  Hence  the  ocean  frequency  response 
H-^  (/)  must  be  computed  by  MMPE  at  this  higher  frequency  (Inkj L+  f^).  The  time- 
sampled  complex  envelope  of  the  received  signal  is  then  given  by 


q[n] 


1  L-l 

Lh  ^ 


^Ink  „  Yr;  1 


.Inkn 


(3.34) 


(Compare  this  to  Eq.  (3.12)).  The  sampled  real  bandpass  signal  at  the  receiver  is  thus 
given  by 
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(3.35) 


3.  Input  Parameters  for  the  MMPE  Model 

The  input  data  for  the  MMPE  model  fall  into  three  main  categories:  the  main 
input  file,  the  source  data,  and  the  environmental  data.  Each  of  these  is  described  below. 

The  main  input  file  is  called  pef  iles  .  inp:  it  is  the  only  filename  that  must 
exist.  All  other  input  files  are  specified  in  this  file.  Eigure  3.6  shows  that  the  first  several 
entries  of  this  main  file  simply  define  the  other  input  files.  Eine  1  specifies  the  source 
data  file.  Eines  2  to  6  specify  files  that  provide  the  environmental  data,  while  line  7 
specifies  the  output  data  file  (a  binary  file).  Eollowing  these  data  files,  the  user  specifies 
the  number  of  points  in  depth,  minimum  depth,  maximum  depth  (at  line  8),  the  number  of 
points  in  range,  minimum  range,  maximum  range  (at  line  9).  The  final  line  in  this  file 
represents  the  vertical  EET  size,  the  range  step,  the  maximum  depth  of  the  calculation, 
and  the  reference  sound  speed. 


Environmental 

data 


D  1  i /► 


-< 

> 


stack:  Base 


Name  of  source  data  input  file 

ipesrc.inp 

i 

Nane  of  ssp  data  input  file 

:pessp.inp 

Name  of  bottoi  bathy  data  file 

:pebath.inp 

Name  of  bottom  properties  file 

:pebotprop.inp 

Name  of  deep  bottom  bathy  file 

:pedbath,inp 

Name  of  deep  bottom  props  file 

ipedbotprop.inp 

Name  of  output  data  file 

:posfc400bw200c.bin 

nzout,  depmin[m],  depmax[m] 

:128  0  200 

nrout,  rnpin[km],  rnpax[km] 

:1000  0  10.0 

nz,  dr[km],  depcalc[m],  cO[m/s] 

:256  0.0  0.0  1500.0 

<J 

_ 1^ 

Eigure  3.6.  pef  iles  .  inp  Eile  of  the  Main  Input  Eile. 


The  source  data  input  file  depicted  in  Eigure  3.7  defines  all  of  the  source 

information.  There  are  two  source  types  which  are  available:  a  wide-angle  source  which 

approximates  a  point  source,  and  a  vertical  line  array  which  approximates  a  continuous 

line  array.  Eor  this  work,  we  will  use  the  wide-angle  point  source  approximation  The 
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first  line  in  Figure  3.7  represents  the  source  cfepth.  The  second  line  denotes  the  array 
length;  an  array  length  of  zero  indicates  a  point  source.  For  a  wide-angle  source,  the  D/E 
angle  is  ignored  since  a  point  source  cannot  be  steered.  The  remaining  lines  represent  the 
center  frequency,  the  bandwidth,  and  the  number  of  the  frequencies  at  which  the  ocean 
response  is  to  be  evaluated  (size  of  the  DFT).  This  number  must  be  a  power  of  two. 


Figure  3.7.  pesrc  .  inp  File  of  the  Source  Data. 

The  environmental  data  consists  of  five  input  files.  The  pessp  .  i  np  file  contains 
the  sound  speed  profile  data  and  is  illustrated  in  Figure  3.8.  As  shown  there,  the  first  line 
contains  two  numbers  indicating  the  number  of  azimuthal  radials  and  the  total  azimuthal 
aperture.  The  second  line  contains  a  single  number  indicating  the  number  of  sound  speed 
profiles.  The  following  line  denotes  the  range  of  the  current  profile  and  the  number  of 
sound  speed  values.  Finally,  the  profile  is  defined  by  pairs  of  the  depth  and  sound  speed. 

The  other  environmental  files  are  similar  in  structure.  Their  detailed  descriptions 
are  contained  in  the  Ocean  Acoustic  Library  web  page 
http ://oalib . s aic . com/PE/mmpe .html  which  is  supported  by  the  U.S.  Office  of  Naval 
Research.  The  pebath .  inp  tile  contains  the  bathymetry  of  the  water/bottom  interface; 
the  pebotprop .  inp  file  contains  the  acoustics  parameters  of  the  medium  just  below 
the  water/bottom  interface;  the  pedbath  .  inp  file  defines  the  “deep”  layer  bathymetry 
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beneath  the  water/sediment  interface;  and  finally,  the  pedbotprop .  inp  file  contains 
the  acoustic  properties  of  the  deep  layer. 


The  last  file,  the  output  data  file  is  a  binary  file  produced  by  MMPE.  It  is  read  by  the 

postprocessing  program  to  produce  final  results 

C.  OCEAN  ENVIRONMENTAL  CHARACTERIZATION 

For  our  simulation  studies  of  underwater  communication,  we  focused  on  shallow 
water  environments.  Three  different  ocean  environments  which  have  three  different 
sound  speed  profiles  (SSP)  were  simulated  to  test  the  performance  of  underwater 
communication  using  a  BPSK  signal.  The  complete  set  of  MMPE  input  files  for  all  three 
cases  are  given  in  Appendix  A. 

I.  Case  I:  Positive  SSP  Gradient 

For  this  case,  the  source  is  located  at  a  range  of  0  m  and  at  a  depth  of  5  m.  The 
sound  speed  is  range- independent  with  a  positive,  linear  gradient  SSP  of  1497  m/s  at  the 
surface  and  1499  m/s  at  a  depth  of  100m.  The  density  of  water  is  assumed  to  be  1.0 
g/cm  .  The  bottom  is  chosen  to  be  a  ‘bumpy’  bottom  with  the  following  depths:  102m  at 
0  km;  101  m  at  2.5  km;  99m  at  5  km;  102  m  at  7.5  km;  and  99  m  at  10  km.  Its 

compressional  sound  speed  is  1700  m/s,  the  sound  speed  gradient  I  density  1.5 
g/crn^,  and  compressional  attenuation  0.1  dB/km/Hz.  Figure  3.9  shows  the  positive 
sound  speed  profile  versus  depth  and  the  corresponding  sound  transmission  loss  plot 
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(versus  range  and  depth)  for  a  source  at  depth  5  m  and  frequency  of  400  Hz. 
Transmission  loss  is  defined  as 


rL  =  201og[p(l,  zj/p(r,z)] 


(3.36) 


where  p{r,z)  is  the  acoustic  pressure  amplitude  and  p(l,zj  is  the  reference  acoustic 
pressure  amplitude  measured  at  range  1  m  (from  the  source)  and  at  the  source  depth  . 

The  color  in  Figure  3.9  represents  the  transmission  loss  in  dB  according  to  the  scale 
shown  on  the  right. 


Figure  3.9.  (a)  Positive  SSP  (b)  Sound  Propagation  at  Frequency  400  Hz,  Source  5  m. 


2.  Case  2:  Strong  Negative  SSP  Gradient 

For  this  case,  the  source  is  located  at  a  range  of  0  m  and  at  a  depth  of  30  m.  The 
sound  speed  is  range -independent  with  a  strong  negative,  bilinear  gradient  (downward 
refraction)  sound  speed  profile  of  1528  m/s  at  the  surface  and  1510  m/s  at  a  depth  of 
50m,  and  1489  m/s  at  the  depth  of  100  m.  The  density  of  water  is  assumed  to  be  1.0 
g/cm  .  The  bottom  is  chosen  to  be  a  ‘bumpy’  bottom  with  the  same  depths  as  in  the 
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previous  case.  Its  compressional  sound  speed  is  1700  m/s,  the  sound  speed  gradient  1 
5“',  density  1.5  g/cm^,  and  compressional  attenuation  0.1  dB/km/Hz.  Figure  3.10  shows 
the  strong  negative  sound  speed  profile  versus  depth  and  the  sound  transmission  loss  for  a 
source  at  depth  30  m,  and  frequency  of  400  Hz. 


tcuKj  |nV9| 


Figure  3. 10. (a)  Negative  SSP  (b)  Sound  Propagation  at  Frequency  400  Hz,  Source  30  m. 


3.  Case  3:  Negative  SSP  Gradient  Below  Surface  Duct 

For  this  case,  the  source  is  again  located  at  a  range  of  0  m  and  at  a  depth  of  5  m. 
The  sound  speed  is  range -independent  with  a  negative,  bi- linear  gradient  sound  speed 
profile  with  surface  duct  of  1492  m/s  at  the  surface  and  1500  m/s  at  a  depth  of  40m,  and 
1489  m/s  at  the  depth  of  100  m.  The  density  of  water  is  assumed  to  be  1.0  g/cm  .  The 
bottom  is  chosen  to  be  a  ‘bumpy’  bottom  with  the  same  depths  as  in  the  previous  cases. 

Its  compressional  sound  speed  is  1700  m/s,  the  sound  speed  gradient  1  5“' ,  density  1.5 
g/cm,  and  compressional  attenuation  0.1  dB/km/Hz.  Figure  3.11  shows  the  positive 
sound  speed  within  the  surface  duct,  the  negative  sound  speed  profile  below  the  surface 
duct,  and  the  sound  transmission  loss  for  a  source  at  depth  5  m  and  frequency  of  400  Hz. 


35 


Figure  3.11.  (a)  SSP  (b)  Sound  Propagation  at  Frequency  400  Hz,  Source  5  m. 

D.  BPSK  DEMODULATION  AND  DETECTION 

At  the  receiver  location,  the  received  signal  waveforms  arrive  distorted  due  to 
noise  and  intersymbol  interference  in  the  underwater  communication  channel.  We  need  to 
demodulate  the  received  signal  in  order  to  recover  the  transmitted  signal  and  recover  the 
binary  data.  Figure  3.12  shows  two  basic  steps  in  the  demodulation  and  detection  process. 
Step  1,  the  waveform- to- sample  transformation,  consists  of  the  demodulator  followed  by 
a  sampler.  At  the  end  of  each  symbol  time  duration  T^ ,  the  sampler  produces  an  output 

z(T^)  which  in  the  absence  of  noise  is  proportional  to  the  energy  of  the  received  symbol. 
Step  2  is  a  decision  process  where  z(Tj,)  is  compared  to  a  threshold  Yo  to  decide  if  the 
received  data  represents  a  binary  1  or  binary  0. 

1.  Correlation  Receiver 

As  mentioned  in  Chapter  II,  there  are  two  main  types  of  degradation  factors  in  the 
performance  of  underwater  communication,  namely  noise  and  dispersion  introduced  by 
the  ocean  channel  with  its  multipath  environment.  In  this  section,  we  ignore  the 
degrading  factor  produced  by  the  ocean  impulse  response  and  assume  that  the  only 
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performance  degradation  is  due  to  AWGN  with  zero  mean.  In  next  section,  the  effect  of 
the  ocean  impulse  response  will  be  added  to  the  received  signal. 


Figure  3.12.  Two  Basic  Steps  in  the  Demodulation/Detection  of  the  Received  Signal. 

In  the  process  of  demodulation,  the  received  signal  is  reduced  to  a  single  random 
variable  z(T^)  sometimes  called  a  “detection  statistic”.  In  the  absence  of  distortion  due 
to  the  ocean,  the  received  signal  can  be  expressed  as: 


q{t)  =  s  (t -x)  +  n{t)  0<t<T 


(3.37) 


where  ^(t)  represents  the  received  signal,  s(t-x)  is  the  transmitted  delayed  BPSK 
signal,  and  n(t)  is  AWGN  with  zero  mean  and  variance,  . 

As  shown  in  Figure  3.13,  the  recovered  signal  is  formed  by  multiplying  the  input 
signal  by  two  local  sinusoidal  carriers  s^it)  and  ^2(0  »  assumed  synchronized  with  the 
received  signal  to  obtain 


y^(t)  =  q(t)-Si(t)=  q(t}  cos(2ti  fj)  (3.38) 
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and 


3^2  (0  =  qit)-S2it)  =  -^(0-  cos(27t  //)  (3.39) 


Figure  3.13.  Correlator  Receiver. 


The  recovered  signals  and  >’2(0  integrated  over  the  bit  interval  to 


obtain 


=  j  =  j  q{t)s.{t)dt  i  =  l,2  (3.40) 


and  the  outputs  of  the  integrations  are  subtracted  to  form  the  detection  statistic 


2(TJ=  Zi(TJ-Z2(TJ 


(3.41) 
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In  the  decision  stage,  z(T^)  is  compared  to  an  optimum  threshold  as  follows 

\  (3.42) 

where  Vj  is  the  signal  component  of  z(T^)  when  s^{t)  is  transmitted,  and  V2is  the  signal 
component  of  z(T^)  when  S2it)is  transmitted.  The  threshold  level  Yq defined  by 
Yo  =  (^1  +  ^2)/3  is  the  optimum  threshold  for  minimizing  the  probability  of  bit  error  [Ref 
1].  For  equal  energy,  equally  likely  antipodal  signals,  where  s^{t)  =  -s^it)  and  v  i=-V2 , 
the  decision  rule  becomes 

z(T,)  Jo=0  (3.43a) 

Assuming  that  s^{t)  corresponds  to  binary  1  and  J2(0  is  binary  0,  the  decision 
rule  thus  reduces  to 

decide  binary  1  if  Zj(T^)  >  Z2(Ti)  (3.43/?) 

decide  binary  0  otherwise 

E.  MULTIPATH  MITIGATION 

In  this  section,  we  assume  the  received  signal  is  further  distorted  by  intersymbol 
interference  due  at  least  in  part  to  multipath  propagation  in  the  channel.  We  need  to  use 
some  multipath  mitigation  technique  to  compensate  for  the  degraded  signal.  One  such 
technique  is  time  reversal.  In  this  method  the  received  signal  is  convolved  with  the  time- 
reversed  impulse  response  of  the  ocean  before  applying  it  to  the  demodulation.  As  shown 
in  Figure  3.14,  the  received  signal  has  the  following  form 

q(t)  =  s(t)*h{t)  +  n{t)  (3.44) 
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where  h{t)  is  the  ocean  impulse  response  in  the  time- domain. 


Local  sinusoidal  carrier 


5j(t)-52(0  =  2cOS(27t/^t) 


Received  signal  Recovered  signal 


A  A 

q(t)  =s{t)*h{t)  +  n{t)  y{t)  -  s{t)- 2cos(2k  fj) ,  where  s{t)  =  s{t)*  h{t)*h  {-t) 


Figure  3.14  Multipath  Mitigation  for  the  Distorted  Signal  in  Underwater  Environment. 

The  distortion  in  the  signal  (represented  by  the  convolution  s[t)*h[t))  is 

mitigated  by  applying  the  time-reversed  ocean  impulse  response  to  the  received  signal.  In 
the  absence  of  noise,  the  resulting  mitigated  signal  is 

s{t)  =  q{t)*h{-t)=  s{t)*h[t)*h  (-t)  (3.45) 

The  convolution  of  h{t)  with  h[-t)  produces  an  impulse- like  signal.  Hence  the 
result  of  Eq.  (3.45)  is  a  tendency  to  restore  the  distorted  signal  to  its  original  condition. 
The  recovered  signal  is  then  formed  by  multiplying  the  mitigated  signal  by  the  two  local 
sinusoidal  carriers  s^it)  and  ^2(0  to  obtain 
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(3.46) 


y(t)=s(l)(s,(t)-s^(t))  =  s(t)-2cos(2ltfJ) 


The  recovered  signals  y(t)  are  integrated  over  the  bit  time  duration  to  form  the 
detection  statistic.  This  is  compared  to  a  threshold  as  before  to  produce  the  recovered 
binary  data. 
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IV.  SIMULATION  RESULTS 


A.  EVALUATION  OF  BIT  ERROR  PROBABILITY  FOR  BPSK  SIGNAL 

In  this  section,  the  parameters  (bandwidth,  sampling  frequency,  bit  rate,  samples 
per  bit,  interpolation  factor  and  up- sampling  frequency)  for  the  BPSK  signal  and  a  Finite 
Impulse  Response  (FIR)  filter  are  defined,  and  the  effect  of  the  AWGN  on  the  bit  error 
performance  is  evaluated. 

I.  Evaluation  for  BPSK  Parameters 

An  important  theorem  of  communication  is  based  on  the  assumption  of  a  strictly 
bandlimited  channel,  i.e.  one  in  which  no  signal  power  whatever  is  allowed  outside  the 
band  of  interest.  For  our  work,  we  need  to  define  the  bandwidth  for  the  BPSK  signal 
transmission.  The  single- sided  power  spectral  density  for  the  BPSK  signal  (also  known  as 
the  analytic  signal)  is  given  by 


‘  BPSK 


if) 


A(T 


siuTt  (/-/,)U 


(4.1) 


This  follows  from  Rj.  (3.8)  by  dropping  the  terms  for  negative  frequencies.  This 

power  spectral  density  is  depicted  in  Figure  4.1  and  is  seen  to  consist  of  a  main  lobe  and 
smaller  sidelobes.  Although  there  are  many  criteria  for  measuring  bandwidth,  for  our 
digital  communication,  we  are  constrained  to  two  bandwidth  criteria  namely  Null-to-null 
bandwidth  and  Power  bandwidth. 

The  Null  -to- null  bandwidth  is  given  by  W  =  2R  =  2/T ^  where  R  is  the 

position  of  the  first  null  relative  to  the  center  frequency  (see  Figure  4.1).  The  sampling 
frequency  corresponding  to  this  definition  of  bandwidth  is  given  by 

f,=2R  =  2IT,  (4.2) 
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Thus,  the  bit  rate  for  this  bandwidth  definition  is  i?  =  iy/2  =  1/T^  and  the  number 
of  samples  per  bit  is  =  2  . 


Figure  4.1.  The  Power  Spectral  Density  for  BPSK  Signal. 


The  Power  bandwidth  defines  the  frequency  band  in  which  99%  of  the  total 
power  resides.  This  bandwidth  has  been  adopted  by  the  Federal  Communication 
Commission  (FCC  Rules  and  Regulations  Section  2.202)  and  states  that  the  occupied 
bandwidth  is  the  band  that  leaves  exactly  0.5%  of  the  signal  power  above  the  upper  band 
limit  and  exactly  0.5%  of  the  signal  power  below  the  lower  band  limit.  Thus  for  this 
definition  99%  of  the  signal  power  is  inside  the  band.  For  the  BPSK  signal  this 
bandwidth  is  given  by  IT  =  20.56R  =  20.56/T^  .  The  sampling  frequency  for  the  power 
bandwidth  is  given  by 


/,  =  20.56R  =  20.56/T, 
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(4.3) 


Thus,  the  bit  rate  for  this  bandwidth  is  R  =  W'/20.56  =1/T^  and  the  number  of 
samples  per  bit  is  =  20.56  . 

Since  the  carrier  frequency  is  generally  much  higher  than  the  bandwidth  of  the 
baseband  signal,  a  random  binary  signal  must  be  interpolated  to  reconstruct  a  bandlimited 
waveform  without  error.  The  interpolation  factor  is  determined  as 


I  =  ceil 


2(/.+iy) 

/, 


(4.4) 


where  I  represents  the  interpolation  factor  and  ceil  rounds  towards  plus  infinity.  The 
increased  sampling  rate  is  then 


(4.5) 

where  represents  the  up- sampling  frequency. 

The  interpolated  random  binary  signal  must  then  be  filtered  to  remove  the 
unwanted  spectral  energy  above  the  band.  A  digital  Finite  Impulse  Response  (FIR)  filter 
is  used.  The  FIR  filter  is  designed  using  a  Hamming  window  with  the  following 
parameters. 


1.0(271  iy/2) 

|(D  „  = - ; - =  passband  digital  frequency 

f s 

2.0(27liy/2)  ^ 

ICO.,  = - ^  =  Stopband  digital  frequency 


\N,= 


87t 


CO  -CO 

s  p 


i  filter  order  (number  of  coefficients  is  + 1) 
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The  filter  parameters  are  shown  for  a  general  FIR  filter  designed  using  the 
Hamming  window  method  in  Figure  4.2.  The  variable  CO  ,  is  the  digital  frequency  defined 
as 


(D  = 


2Tlf 


(4.6) 


Figure  4.2.  General  FIR  Filter  (Hamming  Window  Design). 
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In  the  following  experiments,  the  frequency  response  of  the  BPSK  signal  for  the 
Null- to- null  bandwidth  and  Power  bandwidth  is  shown.  Figure  4.3  shows  the  frequency 
response  for  1000  bits  of  a  BPSK  signal  which  is  sampled  according  to  the  Null-to-null 
bandwidth  and  has  parameters  7?  =100  bits/sec,  /.=400Hz,  =200  Hz.  N^^=2,  1=6 

and  =1200  Hz.  Figure  4.4  shows  the  frequency  response  for  1000  bits  of  a  BPSK 
signal  which  is  sampled  according  to  the  Power  bandwidth  criterion  with  parameters 

7?=100  bits/sec,  /,=4000  Hz,  /,=2100  Hz,  1  =  6,  /,'  =  12600  Hz  and  A^,i=21 
=ce// (20.56) .  The  complete  procedures  to  generate  the  BPSK  signal  are  explained  in 
Appendix  B. 


Figure  4.3.  Frequency  Response  for  BPSK  Signal  Sampled  by  Null-to-null 

Bandwidth. 


As  shown  in  Figure  4.3  and  Figure  4.4,  sampling  according  to  the  Power 
bandwidth  is  more  reliable  than  the  Null-to-null  bandwidth  since  99%  of  the  signal  power 
is  within  the  bandwidth.  However,  since  the  most  energy  is  contained  within  the  Null-to- 
null  bandwidth,  it  is  sufficient  to  work  with  the  Null-to-null  bandwidth  for  our  work.  In 
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the  next  subsection,  the  bit  error  performance  under  AWGN  for  BPSK  signal  having 
Null- to- null  bandwidth  is  evaluated. 

2.  Influence  of  AWGN  on  Bit  Error  Probability 


For  evaluation  of  the  bit  error  performance  under  AWGN,  the  bit  error  probability 
was  tested  for  various  values  of  the  average  signal  power  to  average  noise  power  ratio 
( SNR ).  The  average  signal  power  was  kept  constant. 

E 

The  average  signal  power  to  average  noise  power  ratio  is  defined  as  SNR  = — ^ 

At 


where 


b-  ^  is  ths  average  signal  power,  is  the  white  Gaussian  noise  variance. 


and  At  =  — — .  The  SNR  is  thus  given  by 

N., 


SNR  = 


2a„'At 


(4.7) 


or  in  decibels  as 


^  A  ^ 

SNR  (dB)=101og  (SNR)  =  lOlog  ^  / 

At  ^ 


(4.8) 


The  theoretical  value  of  the  bit  error  probability  Pg  from  Eq.  (2.19)  is: 


Ps  =  Q 


J 


e(^/™)=e 


^  '  a3t.  ^ 


2o„  At 

V  ’  °  y 


(4.9) 


48 


In  the  following  experiment,  the  parameters  for  the  simulation  to  evaluate  the  bit 
error  performance  under  AWGN  for  a  BPSK  signal  are  as  follows:  =10000  bits, 

7?  =100  bits/sec,  IP  =200  Hz,  =400  Hz,  =200  Hz,  =2.  The  noise  variance  is 

varied  from  0.01  to  2  in  increments  of  0.1  and  from  3  to  60  in  increments  of  2. 


BER  vs  SNR  for  BRSK  :  Nbit=10000.  Nsb=2,  fs=200  Hz.  fc=  400  Hz 


Figure  4.5.  BER  versus  SNR  (dB)- 


The  experimental  bit  error  rate  estimates  for  various  values  of  SNR  for  the  Null- 
to-null  bandwidth  is  shown  in  Figure  4.5.  It  can  be  seen  that  the  experimental  values 
follow  the  theoretical  values  given  by  Eq.  (4.9)  but  are  over  all  slightly  lower.  This  is 

due  to  the  effect  of  the  FIR  filter  which  removes  some  of  the  noise.  Since  the  effect  of 
AWGN  on  BER  for  BPSK  signal  is  consistant  with  the  theoretical  values,  our  concern  in 
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ocean  communication  will  be  now  with  the  degradation  due  to  ocean  impulse  response 
and  the  multipath  than  the  degradation  due  to  the  additive  noise. 

B.  BIT  ERROR  DEGADATION  AND  MULTIPATH  MITIGATION  FOR 

BPSK  SIGNAL  IN  A  SHALLOW  WATER  ENVIRONMENT 

In  this  final  section,  the  bit  error  performance  degradation  due  to  the  ocean 
impulse  response  is  evaluated.  This  impulse  response  is  generated  from  the  MMPE 
model  in  the  three  different  ocean  environments  as  described  in  Chapter  III.  In  all  three  of 
these  cases,  we  used  the  same  bandwidth  (200  Hz)  and  center  frequency  (400  Hz)  for  the 
BPSK  signal.  In  all  cases  the  performance  is  compared  for  two  situations.  First,  the 
received  signal  is  applied  directly  to  the  demodulation  with  no  prior  mitigation  steps  and 
the  BER  is  evaluated  as  a  function  of  range  and  depth  from  the  transmitter.  Secondly,  the 
received  signal  is  convolved  with  the  time-reversed  ocean  impulse  response  before 
applying  it  to  the  demodulation.  This  mitigation  step  tends  to  compress  the  spreading 
caused  by  the  ocean  impulse  response.  The  BER  is  then  evaluated  and  compared  to  that 
of  the  unmitigated  situation. 

I.  Bit  Error  Performance  Results  for  a  Positive  SSP  Gradient 

The  parameters  of  this  environment  are  given  in  Chapter  III  (see  Figure  3.9). 
First,  by  choosing  an  ocean  impulse  response  extracted  from  the  MMPE  model  at  a 
chosen  depth  of  50  m  and  range  5  km,  the  influence  of  the  ocean  impulse  response  on  the 
bit  error  rate  at  this  chosen  depth  and  range  is  evaluated.  Then  the  time-reversed  ocean 
impulse  response  is  convolved  with  the  received  signal  to  compensate  for  the  distorted 
signal  and  the  bit  error  rate  is  evaluated  again.  In  the  following,  we  describe  this  one 
specific  case  in  detail.  Eater,  in  section  B  the  evaluation  is  performed  at  many  chosen 
depths  and  ranges  of  this  ocean  environment. 

a.  Investigation  Of  The  One  Specific  Case 

The  parameters  of  the  generated  random  binary  data  are  as  follows: 

=20  bits,  R=100  bits/sec,  Vj^=24,  /^=2400  Hz.  This  discrete  random  binary  data 

b\n\  is  filtered  to  remove  an  unwanted  signal  (i.e.,  above  bandwidth,  W  =200  Hz)  by  a 

FIR  filter.  The  filtered  random  binary  data  is  [n]  =  /?  [n]  [n]  (see  Appendix  B). 

This  filtered  random  binary  data  [  n\  is  modulated  by  cosine  modulating  signal, 
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cos 


^  271  /  n  ^ 

- —  ,  having  a  carrier  frequency  of  /^=400  Hz  and  the  sampling  interval 

V  J 

A7  =  1//^  .  This  entire  procedure  simulates  the  sampling  of  an  analog  BPSK  waveform 
(which  is  what  would  actually  be  present  in  the  water)  at  2400  samples/sec.  The 
modulated  BPSK  signal,  as  shown  in  Figure  4.6  in  time  and  frequency  domains,  is 
271 


>[n]  =  Z7^  [n]  cos 


fs 


The  modified  oeean  frequency  response  H  [f)  extraeted  from  the 

MMPE  model  at  a  ehosen  depth  of  50  m  and  range  5  km  eoriesponds  to  a  eenter 
frequency  of  =  400  Hz,  a  bandwidth  of  W  =300  Hz  and  a  sampling  frequency  of 

/j  =  300  Hz.  To  obtain  the  oeean  impulse  response  at  the  same  sampling  rate  and 
bandwidth  as  the  BPSK  signal,  the  ocean  impulse  response  from  MMPE  is  up- sampled 


by  a  factor  of  8, 


7  = 


2(3x/J 


=  8  ,  and  filtered  to  remove  an  unwanted  signal  (i.e.,  above 


bandwidth,  200  Hz)  by  a  FIR  filter  (see  Appendix  C).  Figure  4.7  shows  the  passband 
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ocean  frequency  response  H^[f)  (magnitude  and  phase)  and  ocean  impulse  response 
[n] .  The  complete  procedure  to  generate  the  passband  ocean  impulse  response  at  this 
high  sampling  frequency  is  explained  in  Appendix  C. 


PafifcUiirld  OCOAn  l-rACi  RAApOriAA 


^  1  ^ 
I  . 


&  3 
Q 


p-U 


r-cqueocy  (Hz) 

©f  Oc©»»>  r»«<<  Ro«f>oA>e 


- 1 

- 1 

- 1 

1 

1 

.  X.. 

— 

. X . 

- 

. 

1 

1 

t-rAgoancy  |Hi) 


Passband  Ocean  Impulse  Response 

T 


Figure  4.7.  Passband  Ocean  Frequency  Response  (Magnitude  and  Phase)  and  Ocean 

Impulse  Response. 


The  BPSK  signal  s[n]  is  convolved  with  the  ocean  impulse  response 
K[^]  by  linear  convolution  in  the  time  domain.  The  received  signal  is 
^[n]  =  5[n]  .  This  simulates  the  distortion  caused  by  propagation  through  the 
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ocean.  Figure  4.8  shows  the  distorted  received  signal  in  time  and  frequency  domains  at 
the  chosen  receiver  location. 


The  recovered  signal  is  obtained  by  multiplying  the  received  signal  q[n\ 


by  the  cosine  modulating  signal,  cos 
2ti  f^n 


f 

\  J  s  y 


Thus,  the  recovered  signal  is 


y[n]  =<7[n]cos 


fs 


.  This  recovered  signal  is  filtered  to  remove  unwanted 


frequency  components  above  the  bandwidth  by  a  FIR  filter  and  integrated  over  the  bit 
time  duration  to  form  the  detection  statistic.  This  is  compared  to  a  threshold  to  produce 
the  recovered  binary  data.  Figure  4.9  shows  that  the  recovered  binary  data  has  a  very  high 
bit  error  rate  (55  %)  due  to  the  ocean  impulse  response  when  compared  to  the  transmitted 
binary  data. 


Figure  4.8.  Received  Signal  in  Time  Domain  and  Frequency  Domain. 
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Figure  4.9.  Comparison  of  Recovered  Binary  Data  with  Transmitted  Binary  Data  for 

the  Received  Signal  (55  %  error  rate). 


To  mitigate  this  high  bit  error  rate  due  to  the  ocean  impulse  response,  the 
time-  reversed  ocean  impulse  response  is  applied  to  the  received  signal.  The  time- 
reversed  ocean  impulse  /i,  \n\  is  obtained  by  reversing  the  passband  ocean  impulse 

response  /?*[«]  in  time  (i.e.  /i,  [n]  = [-n] ).  Figure  4.10  shows  the  passband  time- 

reversed  ocean  frequency  response  (magnitude  and  phase)  and  ocean  impulse  response. 
The  phase  is  the  negative  of  the  phase  in  Figure  4.7. 

To  show  the  mitigated  ocean  impulse  response,  the  time-reversed  ocean 
impulse  response  is  convolved  with  the  ocean  impulse  response.  The  mitigated  ocean 
impulse  response  is  thus  given  by  =\\n\^h^  \n\  =  -n\  ■  Figure  4. 1 1  shows 

that  in  the  absence  of  noise  the  mitigated  ocean  impulse  produces  an  impulse- like  signal. 
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Figure  4. 10. 


Time  Reversed  Ocean  Frequency  Response  (Magnitude  and  Phase)  and 


Ocean  Impulse  Response. 
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Figure  4.11.  Mitigated  Ocean  Impulse  Response  and  Frequency  Response. 


The  mitigated  signal  is  obtained  by  convolving  \  [-n]  with  the  received 
signal  to  obtain  5  [  n]  =  g  [  n]  *  [-n] .  The  recovered  signal  is  obtained  by  multiplying  the 

mitigated  signal  s[n]  by  cosine  to  obtain  y[n]  =5[n]  cos 

same  procedures  to  recover  the  transmitted  binary  data  described  before.  The  recovered 
binary  data  is  compared  to  the  transmitted  binary  data  to  compute  BER.  Figure  4.12 
shows  the  mitigated  signal  in  time  and  frequency  domains  obtained  by  convolving  the 
time-reversed  ocean  impulse  response  with  the  received  signal.  Figure  4.13  shows  that 
the  mitigated  procedure  was  able  to  reduce  the  bit  error  rate  from  0.55  %  to  0  %. 
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Figure  4.12.  Mitigated  Signal  in  Time  and  Frequency  Domain. 


Figure  4.13.  Comparison  of  Recovered  Binary  Data  with  Transmitted  Binary  Data  for 

Mitigated  Signal  Showing  Zero  Error. 
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b.  Results  For  Many  Chosen  Ranges  and  Depths 

In  the  following  experiment,  to  obtain  more  complete  results  of  simulation 


in  this  ocean  environment,  the  simulation  was  performed  for  the  received  signal  and  the 
mitigated  signal  to  evaluate  the  bit  error  performance  for  many  chosen  depths  and  ranges. 
The  parameters  for  this  simulation  are  the  same  as  in  the  previous  section  except  the 
number  of  bits  was  increased  from  20  to  10,000  =10,000).  The  ocean  frequency 

responses  H  [f]  were  extracted  from  the  MMPE  model  at  depths  of  5  m  to  95  m  in 
increments  of  5  m  and  ranges  of  0.5  km  to  9.5  km  in  increments  of  0.5  km  with  a 
sampling  frequency  =  300  Hz,  bandwidth  W  =300  Hz  and  center  frequency  =400 
Hz.  The  experimental  procedures  followed  were  as  described  before. 

As  shown  in  Figure  4.14,  the  recovered  binary  data  set  for  the  received 
signals  at  chosen  depths  and  ranges  as  described  above  has  a  high  bit  error  rate 
(average=0.4946)  due  to  the  distortion  of  the  transmitted  BPSK  signal  due  to  the 
multipath  propagation  resulting  from  the  positive  SSP  gradient.  Figure  4.15  however 
shows  that  the  recovered  binary  data  set  for  the  mitigated  signals  obtained  by  convolving 
the  time-reversed  ocean  impulse  response  with  the  received  signal  has  a  low  bit  error  rate 
(average=0.0451)  compared  to  the  bit  error  rate  of  the  unmitigated  received  signals.  The 
results  of  this  simulation  at  many  chosen  depths  and  ranges  is  consistent  with  the  result  of 
simulation  at  the  former  chosen  depth  of  50  m  and  range  5  km.  We  conclude  therefore 
that  by  using  the  time-reversed  ocean  impulse  response,  we  can  compensate  for  the 
degradation  in  the  bit  error  performance  due  to  the  multipath  propagation  in  this 
simulated  ocean  environment. 
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Figure  4.14.  Bit  Error  Performance  for  Received  Signal  at  Chosen  Depths  5  to  95  m 
and  Ranges  0.5  to  9.5  km  for  a  Positive  SSP  Gradient. 
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Figure  4. 15.  Bit  Error  Performance  for  Mitigated  Signal  at  Chosen  Depths  5  to  95  m 
and  Ranges  0.5  to  9.5  km  for  a  Positive  SSP  Gradient. 


2.  Bit  Error  Performance  Results  for  a  Strong  Negative  SSP  Gradient 

In  the  following  experiment,  for  a  shallow  water  acoustic  channel  of  having  a 
strong  negative  SSP  gradient,  the  simulation  was  performed  for  the  received  signal  and 
the  mitigated  signal  to  evaluate  the  bit  error  performance  for  many  chosen  depths  and 
ranges.  The  SSP  and  TL  plots  are  shown  in  Figure  3.10a  and  3.10b.  The  signal 
parameters  and  grid  of  range  and  depth  values  were  the  same  as  for  the  previous  case  and 
the  experimental  procedures  followed  were  as  described  before.  As  shown  in  Figure  4.16, 
the  recovered  binary  data  set  for  the  received  signals  at  chosen  depths  and  ranges  as 
described  above  has  a  high  bit  error  rate  (average=0.4981)  due  to  the  distortion  of  the  ISI 
resulting  from  the  strong  negative  sound  speed  profile.  Figure  4.17  shows  that  the 
recovered  binary  data  set  for  the  mitigated  signals  however  has  a  low  bit  error  rate 
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(average=0.0415)  compared  to  the  bit  error  rate  of  the  unmitigated  received  signals.  We 
conclude  therefore  that  by  using  the  time-reversed  ocean  impulse  response,  we  can 
compensate  for  the  degradation  in  the  bit  error  performance  due  to  the  multipath 
propagation  in  this  ocean  environment  as  well. 


Figure  4.16.  Bit  Error  Performance  for  Received  Signal  at  Chosen  Depths  5  to  95  m 
and  Ranges  0.5  to  9.5  km  for  a  Strong  Negative  SSP  Gradient. 
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Figure  4. 17.  Bit  Error  Performance  for  Mitigated  Signal  at  Chosen  Depths  5  to  95  m 
and  Ranges  0.5  to  9.5  km  for  a  Strong  Negative  SSP  Gradient. 

3.  Bit  Error  Performance  Results  for  a  Negative  SSP  Gradient  below 
Surface  Duct 

In  the  following  experiment,  for  a  shallow  water  acoustic  channel  having  a 
negative  SSP  gradient  below  surface  duct,  the  simulation  was  performed  for  the  received 
signal  and  the  mitigated  signal  to  evaluate  the  bit  error  performance  for  many  chosen 
depths  and  ranges.  The  SSP  and  TL  plots  for  this  environment  are  shown  in  Figure  3.11a 
and  3.11b.  The  signal  parameters  and  grid  of  range  and  depth  values  were  the  same  as  in 
the  previous  case  and  the  experimental  procedures  followed  were  as  described  before.  As 
shown  in  Figure  4.18,  the  recovered  binary  without  mitigation  has  a  high  bit  error  rate 
(average=0.4989)  due  to  the  multipath  and  the  resulting  ISI.  Figure  4.19  shows  that  the 
recovered  binary  data  set  for  the  mitigated  signals  however  has  a  low  bit  error  rate 
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(average=0.0398)  compared  to  the  bit  error  rate  of  the  unmitigated  received  signals.  We 
conclude  that  for  this  environment  as  well,  by  using  the  time-reversed  ocean  impulse 
response,  we  can  compensate  for  the  degradation  in  the  bit  error  performance  due  to  the 
multipath  propagation. 
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Figure  4.18.  Bit  Error  Performance  for  Received  Signal  at  Chosen  Depths  5  to  95  m 
and  Ranges  0.5  to  9.5  km  for  a  Negative  SSP  Gradient  below  Surface  Duct. 
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Figure  4. 19.  Bit  Error  Performance  for  Mitigated  Signal  at  Chosen  Depths  5  to  95  m 

and  Ranges  0.5  to  9.5  km  for  a  Negative  SSP  Gradient  below  Surface  Duct. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  primary  goal  of  this  thesis  was  to  investigate  mitigating  the  degradation  on 
the  bit  error  performance  of  the  BPSK  signal  by  convolving  the  time-reversed  ocean 
impulse  response  with  the  received  signal  distorted  by  the  multipath  propagation  in  three 
different  shallow  water  acoustic  channels.  The  simulation  was  validated  by  a  comparison 
of  the  results  of  the  bit  error  performance  from  the  received  signal  and  those  from  the 
mitigated  signal.  In  our  experiments,  it  was  possible  to  reduce  the  high  bit  error  rate  as 
follows.  For  a  positive  SSP  gradient,  the  average  bit  error  rate  decreased  from  0.4946  to 
0.0451.  For  a  strong  negative  SSP  gradient,  the  average  bit  error  rate  decreased  from 
0.4981  to  0.0415.  For  a  negative  SSP  gradient  below  the  surface  duct,  the  average  bit 
error  rate  decreased  from  0.4989  to  0.0398.  Thus,  it  was  possible  to  improve  the  bit  error 
performance  of  a  BPSK  signal  by  using  the  time-reversed  ocean  impulse  response. 

The  experiments  assumed  that  the  correct  ocean  impulse  response  was  used  at 
each  position  of  the  receiver.  The  sensitivity  to  range  and  depth  or  incorrect  ocean 
impulse  response  was  not  investigated  in  any  quantitative  manner.  There  were  some 
indications,  however,  of  a  fair  amount  of  sensitivity  to  change  in  position. 

B.  RECOMMENDATIONS 

Although  our  experiment  was  not  able  to  produce  a  sufficiently  low  bit  error  rate 
needed  to  achieve  an  effective  underwater  communication  in  the  ocean  environment,  the 
algorithm  may  be  of  some  use  if  appropriate  error  correction  coding  is  employed  to 
reduce  the  bit  error  rate.  Experiments  with  a  time-reversed  ocean  impulse  response 
showed  reasonable  success  in  reducing  the  bit  error  rate  of  the  distorted  BPSK  signal  due 
to  the  intersymbol  interference.  The  combination  of  appropriate  other  filtering  and 
coding/decoding  with  this  time-reversed  ocean  impulse  response  may  therefore  be  worthy 
of  further  investigation. 
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APPENDIX  A.  MMPE  INPUT  EILES  EOR  THREE  DIEEERENT 
OCEAN  ENVIRONMENTAL  CASES 


This  appendix  gives  the  complete  set  of  MMPE  input  files  for  three  different 
ocean  environmental  cases  as  described  in  Chapter  III.  The  MMPE  input  files  of  the 
environmental  data  except  the  input  file  of  the  sound  speed  profile  and  pesrc  .  inp  file 
of  the  source  data  are  same  for  three  cases. 

A.  MMPE  INPUT  FILES  FOR  POSITIVE  SSP  GRADIENT 
1.  pef  lies  .  inp  File  of  the  Main  Input  File 


□  1^  H  #  '  II  ^  M  /►  1 5^  ^  ir  ii  jii 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Name  of  source  data  input  file 
Name  of  ssp  data  input  file 
Name  of  bottom  bathy  data  file 
Name  of  bottom  properties  file 
Name  of  deep  bottom  bathy  file 
Name  of  deep  bottom  props  file 
Name  of  output  data  file 
nzout,  depmin[m],  depmax[m] 
nr out,  rngmin[km],  rngmax[km] 
nz,  dr[km],  depcalc[m],  c0[m/3] 


: pesrc. inp 
:pessp. inp 
:pebath. inp 
ipebotprop. inp 
ipedbath. inp 
:pedbotprop. inp 
: posf c400bw300b . bin 
:10  0  200 
:10  0  10.0 
:256  0.0  0.0  1500.0 


Stack:  jBase 


i\ 


2.  pesrc .  inp  File  of  the  Source  Data 
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3 


pessp .  inp  File  of  the  Environmental  Data 


D 

S|  ^ 

1 

1  0. 

2 

1 

3 

0.  101 

4 

0 

1497 

5 

1 

1497 

6 

2 

1497 

7 

3 

1497.1 

8 

4 

1497.1 

9 

5 

1497.1 

10 

6 

1497.1 

11 

7 

1497.1 

12 

8 

1497.2 

13 

9 

1497.2 

14 

10 

1497.2 

15 

11 

1497.2 

16 

12 

1497.2 

17 

13 

1497.3 

18 

14 

1497.3 

19 

15 

1497.3 

20 

16 

1497.3 

21 

17 

1497.3 

22 

18 

1497.4 

23 

19 

1497.4 

24 

20 

1497.4 

25 

21 

1497.4 

A 

o  rx 


/►I  ®  I 


stack:  Base  .2S] 


4.  pebath .  inp  File  of  the  Environmental  Data 
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5. 


pebotprop .  inp  File  of  the  Environmental  Data 


6.  pedbath .  inp  File  of  the  Environmental  Data 


7.  pef  lies  .  inp  File  of  the  Environmental  Data 
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B.  MMPE  INPUT  FILES  FOR  STRONG  NEGATIVE  SSP  GRADIENT 
I.  pef  lies  .  inp  File  of  the  Main  Input  File 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Name  of  source  data  input  file 
Name  of  ssp  data  input  file 
Name  of  bottom  bathy  data  file 
Name  of  bottom  properties  file 
Name  of  deep  bottom  bathy  file 
Name  of  deep  bottom  props  file 
Name  of  output  data  file 
nzout,  depmin[m],  depmax[m] 
nr out,  rngmin[km],  rngmax[km] 
nz,  dr [km],  depcalc[m],  cO[m/s] 

_ 


:pesrc. inp 
:pessp. inp 
:pebath. inp 
rpebotprop. inp 
:pedbath. inp 
ipedbotprop.  inp 
: negf c400bw300a.  bin 
:128  0  200 
:1000  0  10.0 
:256  0.0  0.0  1500.0 


Stack:  jBase  .2S\ 


JaT 


2.  pessp .  inp  File  of  the  Environmental  Data 


D  0^  a 


/►I  il  I  ^ 


stack:  Base  .2S] 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


<1 


0. 

101 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


0  1528 

1527.6 

1527.3 
1526.9 

1526.6 

1526.2 

1525.8 

1525.5 

1525.1 

1524.8 

1524.4 
1524 

1523.7 

1523.3 
1523 

1522.6 

1522.2 

1521.9 

1521.5 

1521.2 

1520.8 

1520.4 
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C.  MMPE  INPUT  FILES  FOR  NEGATIVE  SSP  GRADIENT  BELOW 
SURFACE  DUCT 

1.  pef  lies  .  inp  File  of  the  Main  Input  File 


D  C?  H  #  #4  /►  I  §]  jgl  I  <s|ll  nr  lil  Jgl  I  stack: 


1 

Name  of  source  data  input  file 

rpesrc. inp 

A. 

2 

Name  of  ssp  data  input  file 

rpessp. inp 

3 

Name  of  bottom  bathy  data  file 

:pebath. inp 

4 

Name  of  bottom  properties  file 

rpebotprop. inp 

5 

Name  of  deep  bottom  bathy  file 

rpedbath. inp 

6 

Name  of  deep  bottom  props  file 

:pedbotprop. inp 

7 

Name  of  output  data  file 

: surf c400bw300a. bin 

8 

nzout,  depmin[m],  depmax[m] 

:128  0  200 

9 

nr out,  rngmin[km],  rngmax[km] 

:1000  0  10.0 

10 

11 

nz,  dr [km],  depcalc[m],  cO[m/s] 

:256  0.0  0.0  1500.0 

▼ 

^1 

LI 

2.  pessp .  inp  File  of  the  Environmental  Data 


D  0^  y  #  1 

^  /► 

t  1 

stack:  j 

jBase  2< 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


1  0. 

▲ 

1 

0.  101 

0  1528 

1 

1527.6 

2 

1527.3 

3 

1526.9 

— 

4 

1526.6 

5 

1526.2 

6 

1525.8 

7 

1525.5 

8 

1525.1 

9 

1524.8 

10 

1524.4 

11 

1524 

12 

1523.7 

13 

1523.3 

14 

1523 

15 

1522.6 

16 

1522.2 

17 

1521.9 

18 

1521.5 

19 

1521.2 

20 

1520.8 

21 

1520.4 

▼ 

4l  I  kl 
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APPENDIX  B.  DETAILED  SIGNAL  PROCESSING  STEPS 


This  appendix  gives  the  complete  procedures  from  the  generation  of  the  BPSK 
signal  to  the  demodulation  of  the  received  signal  for  an  ideal  digital  communication 
channel  (i.e.,  without  any  distortion) 

A.  RANDOM  BINARY  DATA  GENERATION 

The  generated  random  binary  data  b\n\  has  parameters  =  20  bits,  R  =100 

bit/sec,  N =  24 ,  f^  =  2400  Hz  and  T^  =  0.01  sec.  Figure  B- 1  shows  the  random  binary 
data  in  time  and  frequency  domains. 


Birsr,'  C'3t»-  tJbil?  =35 .  tJsb  =34.  ■b=  2430  Hi 


Tin«  (ssc) 

Frrq  Response  of  Binarv  Dola 


Figure  B- 1 .  Random  Binary  Data  in  Time  Domain  and  Frequency  Domain. 
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B.  FILTERING 


This  random  binary  data  is  filtered  to  remove  an  unwanted  spectral  energy  above 
bandwidth  {W  =  200  Hz).  The  FIR  filter  \n\  is  designed  using  a  Hamming  window 

as  described  in  chapter  4.  The  filtered  binary  data  is  bf\n\  =  b\n\^h^pp  \n\  .  Figure  B-2 
shows  the  frequency  response  of  FIR  filter  h^pp  \n\  and  the  filtered  binary  data  bj\n\ . 


Figure  B-2.  Frequency  Response  of  the  FIR  Filter  and  the  Filtered  Binary  Data. 


C.  GENERATION  OF  BPSK  SIGNAL 


The  filtered  binary  data  b^  [  n\  is  modulated  by  a  cosine  modulating  signal, 


k[n]  =  cos 


2Jlf,n 

fs 


having  a  carrier  frequency  of  =  400  Hz.  The  modulated  BPSK 


signal  is  =  [n]'^[n].  Figure  B-3  shows  the  modulated  BPSK  signal  in  the  time 

and  frequency  domains. 


74 


Figure  B-3.  BPSK  Signal  ^[n]  in  Time  Domain  and  Frequeney  Domain. 


D.  DEMODULATION  OF  BPSK  SIGNAL,  FILTERING 

The  transmitted  BPSK  signal  ^[n]  is  demodulated  by  multiplying  by  a  eosine 
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demodulating  signal,  d  [n]  =  cos 


The  demodulated  signal  is  given  by 


y[n]  =  j  [n]  and  produces  a  double  frequency  component  centered  at  ±800  Hz  as 
shown  in  Figure  B-4.  The  signal  y\n\  is  filtered  again  to  remove  unwanted  spectral 
energy  outside  bandwidth  by  using  same  FIR  filter  described  above.  The  demodulated 
and  filtered  signal  is  [n]  =  y\n\ *  [n] .  Figure  B-4  shows  the  frequency  response  of 

the  demodulated  BPSK  signal  y\n\  before  lowpass  filtering.  Figure  B-5  shows  the 
demodulated  and  filtered  signal  y^  \n\  in  time  and  frequency  domains. 
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Figure  B-4.  Demodulated  BPSK  Signal  y\n\  in  Frequeney  Domain. 


Demodul’il'-d  Lcv*p»»9i?d  BPSK  Siun^l 


PraquancY  Paspori'ia  oi  Dcmodutatad  and  Lcrv^assad  BPSK  Signal 


Figure  B-5.  Demodulated  and  Filtered  BPSK  Signal  in  Time  and  Frequeney  Domain. 


E.  BER  COUNTING 

This  demodulated  and  filtered  signal  [n]  is  integrated  over  the  bit  time 
duration  to  form  the  detection  statistic  and  then  compared  to  an  optimum  threshold  as 
described  in  Eq.  (3.43a)  of  Chapter  III  and  by  this  decision  rule,  we  decide  binary  1  or  0. 
Thus,  the  recovered  binary  data  is  compared  to  the  transmitted  binary  data.  Figure  B-6 
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shows  that  the  bit  error  rate  of  the  recovered  binary  data  for  an  ideal  digital 
communication  signal  goes  to  zero. 


RerrMaffrtl  nifiny  Data  Set  BFR  =  D 


TianamtlCril  Btrt»y  DaiaSwl 


1 - 1 - 

*  1 

- 1 

r - ‘ 

*  1 

— 

1 

_ 

0  ;  4  G  C 

i  10  ir  14  1G  1 

5  ~  36 

Figure  B-6.  Comparison  between  Recovered  Binary  Data  and  Transmitted  Binary 

Data. 
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APPENDIX  C.  COMPLETE  PROCEDURES  TO  GENERATE  THE 
PASSBAND  OCEAN  IMPULSE  RESPONSE  AND  EREQUENCY 
RESPONSE  EROM  THE  MMPE  MODEL 


This  appendix  gives  the  complete  procedures  to  generate  the  passband  ocean 
impulse  response  and  frequency  response  at  a  high  sampling  frequency  from  the  modified 

ocean  frequency  response  H  (f)  extracted  from  MMPE  model. 

A.  MODIFIED  OCEAN  FREQUENCY  RESPONSE  AND  OCEAN  IMPULSE 
RESPONSE  FROM  MMPE  MODEL 

The  modified  ocean  frequency  response  H  [f)  extracted  from  MMPE  model  as 

mentioned  in  Eq.  (3.3 1)  of  Chapter  III  has  IT  =  300  Hz,  /^  =  400  Hz,  =300  Hz 

and  n^.  =256,  where  represents  the  number  of  frequency  components  or  EET  size. 

Eet  us  consider  a  chosen  depth  of  50  m  and  range  5  km  in  the  ocean  environment  having 
a  positive  SSP  gradient.  Eigure  C-1  shows  the  modified  ocean  frequency  response 

H  (f)  and  ocean  impulse  response  h[n]  obtained  by  taking  the  inverse  DET  of 


H'if). 

B.  BASEBAND  OCEAN  IMPULSE  RESPONSE  BY  PADDING  ZEROS 

To  obtain  the  ocean  impulse  response  at  the  same  sampling  frequency  as  the 
BPSK  signal,  we  need  to  pad  zeros  between  lower  part  and  upper  part  of  the  modified 

ocean  frequency  response  H  (/).  In  this  case,  to  obtain  8  times  sampling  frequency 

^  =8/^^  ^  )  of  the  original  sampling  frequency  of  =300  Hz,  we  need  to  pad 


W  W 

zeros  between  -3/^  and  except  for  the  interval  These  zeros  should 


be  equally  spaced  as  A/ 


IT 


Uf-l 


T.1765  Hz.  Eigure  C- 2  shows  the  baseband  ocean 


frequency  response  H^[f)  and  ocean  impulse  response  [n] . 
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Figure  C- 1 . 


Modified  Ocean  Frequency  Response  H  [f]  and  Impulse  Response 


Figure  C-2. 


Baseband  Ocean  Frequency  Response  H^[f)  and  Impulse  Response 
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C.  FILTERING  AND  PRE-ENVELOPE  OF  THE  OCEAN  IMPULSE 
RESPONSE 

Since  the  baseband  ocean  response  W  =  300  Hz,  we  need  to  filter  it  to  have  the 
same  bandwidth  as  the  BPSK  signal  (200  Hz).  The  FIR  filter  h^pp \n\  is  designed  using  a 

Hamming  window  as  described  in  chapter  4.  Figure  C-3  shows  the  frequency  response  of 
the  FIR  filter  and  the  filtered  baseband  ocean  frequency  response. 

The  filtered  baseband  ocean  impulse  response  ^[n\  is  multiplied  by 
exp(-y27l to  obtain  the  pre-envelope  of  the  passband  ocean  impulse  response, 
/i^[n]  = /i/^[n]  exp(-y27l/,t) .  Figure  C-4  shows  the  frequency  response  of  the  pre¬ 
envelope  of  the  ocean  impulse  response  H  ^[f) . 


Figure  C-3.  Frequency  Response  of  FIR  Filter  and  Filtered  Baseband  Ocean 

Frequency  Response. 
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Figure  C-4.  Pre-Envelope  Ocean  Frequency  Response 


D.  PASSBAND  OCEAN  IMPULSE  RESPONSE 

Since  the  pre-envelope  of  the  passband  ocean  impulse  response  is 

[n]  -  \  \n\  +  7^;,  [n]  where  \  \n\  is  the  passband  ocean  impulse  response  and  \n\  is 
the  Hilbert  transform  of  \  \n\  (see  Ref  13).  The  passband  ocean  impulse  response  is 
obtained  by  taking  the  real  part  of  hp\n\  and  by  scaling  by  a  factor  of  2.  The  passband 
ocean  impulse  response  is  defined  as  =  2Re  [n]) .  The  passband  frequency 

response  H^{f)  is  obtained  by  doing  the  Fourier  transform  of  /?;,[«].  Figure  C-5  shows 
the  passband  ocean  frequency  response. 
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